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Abstract 

Optical gas detection in microsystems is limited by the short micron scale optical 
path length available. Recently, the concept of slow-light enhanced absorption has 
been proposed as a route to compensate for the short path length in miniaturized 
absorption cells. We extend the previous perturbation theory to the case of a Bragg 
stack infiltrated by a spectrally strongly dispersive gas with a narrow and distinct 
absorption peak. We show that considerable signal enhancement is possible. As an 
example, we consider a Bragg stack consisting of PMMA infiltrated by 02- Here, 
the required optical path length for visible to near-infrared detection (~ 760 nm) 
can be reduced by at least a factor of 10^, making a path length of 1 mm feasible. 
By using this technique, optical gas detection can potentially be made possible in 
microsystems. 
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1 Introduction 



The integration of optics and microfiuidics on lab-on-a-chip microsystems has 
recently been the topic of much research [1,2,3], partly motivated by its diverse 
applications in chemical and biochemical analysis [4]. The miniaturization of 
chemical analysis systems, however, presents many optical challenges since 
light matter interactions suffer from the reduced optical path length L in lab- 
on-a-chip systems compared to their macroscopic counterparts as illustrated 
in Fig. 1. Mogensen et al. [5] demonstrated that for Beer-Lambert absorption 
measurements in a lab-on-a-chip system, a typical size reduction by two orders 
of magnitude severely reduces the optical sensitivity in an inversely propor- 
tional manner. This drawback is even more pronounced for gas detection in 
microsystems as most gases have a very weak absorption line. Thus, microsys- 
tems may not seem the obvious solution for gas sensing and detection. On 
the other hand, the major reason for being interested in pursuing microsys- 
tem opportunities is of course that many applications call for compact gas 
sensors because of limited sensor space. On the more practical level, we em- 
phasize that lab-on-a-chip systems can both address very minute volumes (no 
flow), but also allows for measurements of larger volumes flowing through the 
microchannels where detection takes place. In that way, spatial variations in 
the gas properties can potentially be mapped onto a time axis and monitored 
through time-traces of the optical signal. This space-time mapping is a general 
approach (employed for a variety of other measurements) which is only pos- 
sible in microchannels supporting laminar flow. For the optical sensitivity, it 
has been proposed, that by introducing a porous and strongly dispersive mate- 
rial, such as a photonic crystal, into the lab-on-a-chip microsystems one could 
potentially solve these problems through slow-light enhanced light-matter in- 
teractions [6,7,8,9]. First experimental evidence of this effect for gas sensors in 
the mid-infrared range was reported by Lambrecht et al. [10]. Here, we extend 
the perturbative theory in Refs. [6,7,8,9] to the case of a dispersive Bragg 
stack infiltrated by a spectrally strongly dispersive gas with a narrow and dis- 
tinct absorption peak. We emphasize the example of oxygen monitoring and 
sensing based on the two distinctive bands near the visible to near-infrared 
light range of 760 nm referred to as the O2 A band [11]. Spectral data are 
available from the HITRAN database [12]. Though, the A band has a very 
weak absorption feature, it offers the potential to establish an optical in situ 
O2 detection for application in many fields including combustion processes [13] 
and fire research [14]. 
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Fig. 1. (a) Schematic of a typical macroscopic gas detection setup with a light-matter 
interaction path length of 10 cm. (b) Schematic of a microscopic gas detection setup 
with a light-matter interaction path length of 100 fim made possible by the presence 
of a Bragg stack inside the microchannel. (c) Liquid/gas infiltrated Bragg stack 
composed of alternating solid and gas layers of thickness bg and bg, respectively. 

2 Beer Lambert absorption 



In gas detection experiments, one typically uses Beer-Lambert absorption to 
determine the concentration of the substance of interest. As the light traverses 
the test chamber the light intensity is exponentially attenuated according to 
the relation 

/ = /oexp(-aL), (1) 

where a is the absorption coefficient, typically a linear function of the concen- 
tration C of absorbing molecules (not to be confused with c; the speed of light 
in vacuum to be introduced at a later stage). Any reduction in the optical 
path length L will penalize the sensitivity of the optical measurement signif- 
icantly. However, by using a porous photonic crystal, it has been suggested 
that one may effectively enhance a by slowing down the light via a photonic 
crystal thereby overcoming the short path length obstacle in microsystems [6]. 
To quantify the enhancement, we shall use the ratio 
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where ape is the absorption coefficient inside the photonic crystal and Og 
the corresponding homogeneous space absorption coefficient associated with 
the gas itself. The subscripts are introduced in order to carefully distinguish 
between absorption taking place in the photonic crystal material and the gas, 
respectively. In this work we consider the particular simple one-dimensional 
example of a Bragg stack and extend previous work [15] by using a more 
realistic model for the gas absorption profile. 
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3 The Bragg stack 



The periodic Bragg stack is iUustrated in panel (c) of Fig. 1 and is composed 
of ahernating layers of gas and solid of thickness bg and feg, respectively. The 
lattice constant is A = bg + bs and for the dielectric functions we allow for 
complex-valued bulk parameters, i.e. eg = e'g + ie'g and = e'^ + ie". For later 
convenience we rewrite the dielectric functions in terms of the bulk refractive 
index n and the damping coefficient a, 



eg = nl 



+ ingag/k, (3a) 



es = nl + iusas/k, (3b) 

where k = 2tt/ X = u/c is the free-space wave vector with A being the corre- 
sponding free-space wavelength and u! the angular frequency. The dispersion 
relation for the Bragg stack is governed by (see Ref. [15] and references therein) 

cos (kA) = F{k), (4a) 

where k is the Bloch wave vector and 

F{k) = cos [yfTgkbg^ cos {^/elkb,) - ^-y= sin {y^gkb^ sin (y/e^kbs) . (4b) 



In the following we explicitly write the Bloch wave vector as k = n' + in", 
where k' and k" are real and imaginary parts, respectively. The imaginary 
part causes an exponential damping of the Bloch states with a corresponding 
attenuation coefficient a^c = 2fi;". Likewise, in the homogeneous case we have 
ag = ke''Lfe' and in this way Eq. (2) now becomes 



7 = X ^ ^ . (5) 
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If 7 is greater than unity, it means that we have enhanced absorption due to 
the geometry of the Bragg stack. The absorption enhancement comes about 
because the group velocity of the light duj/dhi' is reduced by the geometry 
of the Bragg stack [6]. In Eq. (5), the slow light phenomenon manifests itself 
through an increasing imaginary part of the Bloch vector k" as the band 
edge is approached. In deriving Eq. (5), we have separated the two different 
contributions to the absorption from the gas and the solid, thus implicitly 
assuming a small change in attenuation caused by the presence of the gas. 
However, we emphasize that the effect of material dispersion is included in 
the dispersion relation in an exact analytical and non-perturbative way and 
thus our results fulfill the Kramers-Kronig relation by taking the full interplay 
between waveguide and material dispersion into account. 
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Fig. 2. Band diagram for a Bragg stack with eg = 1.0^, = 1.5^, bg = 0.8A, and 
bs = 0.2A, corresponding to a polymer infiltrated by gas. 

4 Examples 



In the following we illustrate the prospects for slow-light enhanced absorption 
by means of examples. 



4-1 Non-absorbing materials 



First, we consider the problem of two non-absorbing materials, say a gas and 
a polymer with eg = 1.0^ and = 1.5^. The corresponding band diagram, 
obtained from a numerical solution of Eq. (4a), is shown in Fig. 2. Two band 
gaps, indicated by yellow shading, are visible near uA/c = 2.8 and 5.5. As can 
be seen, the group velocity approaches zero near the band gap edges while 
the imaginary part of k attains a finite value only inside the band gaps where 
propagation of electromagnetic radiation is prohibited due to the exponential 
damping. 
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Fig. 3. Band diagram for the same problem as in Fig. 2, but with UgK = a^A = 0.1, 
corresponding to a finite bulk damping in the gas and the polymeer, respectively. 

4-2 Frequency-independent weakly absorbing materials 

Next, we consider two weakly absorbing materials, say again the same gas and 
polymer, but now with imaginary contributions to the dielectric functions 
corresponding to damping parameters a^A = a^A = 0.1. As seen from the 
corresponding band diagram in Fig. 3, the main cause is to smear out the 
spectrally pronounced dispersion features. First, it increases n" from zero to 
a finite value outside the band gaps and second it smears out the band gap, 
setting a lower limit on the magnitude of the group velocity. 



4-3 Strongly frequency- dependent absorbing gases 

Having treated the simple case of spectrally uniform absorption, we shall now 
consider a slightly more realistic model for the frequency-dependent absorp- 
tion in gasses. Due to causality and the Kramers-Kronig relation, the presence 
of frequency-dependent absorption is inevitably accompanied by material dis- 
persion in the gas which will play in concert with the waveguide dispersion 
contributed by the geometry of the Bragg stack. To this end, we will now 
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Fig. 4. Band diagram for a Bragg stack with OgA = 0.1, = 1.5, bg = 0.8A, 
bg = 0.2A and €g given by equations 6a-6b with loqA/c = 3.1158, aA/c = 10~^ and 
Ac^/A'^ = W~^. The black solid line (b) shows the results for the gas-infiltrated 
Bragg stack, while the corresponding dashed line (a) is in the absence of gas infil- 
tration. In the k" plot, the red solid line (c) shows the absorption in the absence of 
the Bragg stack. Only the spectral region just above the band gap is emphasized in 
the plots. 

consider a Lorentzian absorption line of the form [16] 

= ^71 2 2 ' (6a) 

where the corresponding real part of the dielectric function follows from the 
Kramers-Kronig relation and is given by 

4H = i + ^ , , "IX' , (6b) 

In this model, the strength A scales with the concentration of molecules while 
a quantifies the spectral line width around the resonance frequency ujq. As an 
example the band diagram for OgA = 0.1, ujqA/c = 3.1158, aA/c = 10~^, and 
Ac^/A"^ = 10"^ is shown in Fig. 4. Note the zoom-in on the first band gap 
and the different scale compared to Figs. 2 and 3. The horizontal dashed black 
line indicates the above-band-edge position of the absorption line while the 
yellow shading indicates the band gap for the absorption- free structure. The 
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red line shows bulk-absorption of the gas itself, thus indicating an absorption 
enhancement of the order 7 ~ 5.5. We emphasize that this value depends 
strongly on ag, a^, bg, and bs and as we shall see in the next section, even 
stronger absorption enhancement is possible by carefully tuning the Bragg 
stack parameters to fit the absorption line of interest. 



5 Application to oxygen detection 

Oxygen molecules, O2, exhibit a distinctive absorption band in the visible 
to near-infrared light range at 760 nm, which is called the Oxygen A band 
or the atmospheric transition of oxygen. These transitions are electric dipole 
forbidden, so only magnetic transitions are allowed, that require a change in 
the rotational quantum number of dN = ±1. Therefore the A-band splits 
in to the R-branch [dN = —1), starting at 759.8 nm with a bandwidth of 
2.2 nm and the P-branch [dN = +1) starting at 762.3 nm with a bandwidth 
of 5.7 nm [17,11]. Though, the A band has a very weak absorption feature, it 
offers the potential to establish an optical in situ O2 detection for a variety of 
applications, ranging from combustion processes [13] and fire research [14] to 
biomedical and environmental systems [18]. For this purpose different spectro- 
scopic methods are used to enlarge the sensitivity of absorption measurements. 
Folding the laser beam several times in a multi reflection cell or using the long 
photon life times of high finesse cavities (cavity ring-down spectroscopy) offer 
an enhancement in sensitivity in respect to an enlargement of the absorption 
path length. In contrast, the method of slow light enhancement is related to an 
increased light matter interaction (though it may be phenomenologically ex- 
plained as an enlarged effective path length), leading to higher absorption cross 
sections. This effect offers the potential for micrometer scaled probe volumes 
which is an important feature for e.g. lab on a chip applications. Typically, 
the absorption of oxygen in ambient air at the 760 nm line is of the order 2.6% 
with an optical path length of around 1 m [13]. At this wavelength, the O2 
absorption spectrum is minimally sensitive towards temperature changes and 
it is fairly influenced by the pressure variation [17,11,13]. The lowest oxygen 
concentration that can be measured is about 20-50 ppm which corresponds 
to transmission change of about 0.0005 [11]. The smallest experimentally de- 
tected absorption is about 3 times larger than the shot noise limit and is 
about 10~^. This corresponds to a concentration of about 8 ppm at 1 m path 
length [19,20]. 

By extracting the spectroscopic data from the published works [13,14,21,22,23], 
we correlated the O2 absorption coefficient, a [see Eq. (1)] as a function of O2 
concentration C in the radiated volume. To lowest order in the gas concentra- 
tion we expect that a(C) = r]C and the linearity of the correlation obtained is 
shown in Fig. 5. The slope of the relatively linear line shows that it is possible 
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Fig. 5. O2 absorption coefficient q as a function of the O2 concentration C in the 
radiated volume. The inset shows the A band based on HITRAN data [12] for a 
path length of 1 m. 

to obtain a fairly good measurement of O2 molecule through the excitation 
of light at A ~ 760 nm. From a linear fit we find that 77 ~ 3.7 (mol/L)~^m~^. 
We note that the different absorption lines in the A-band have different line 
strengths which may cause the scatter of data in Fig. 5. The absorption in 
PMMA, at frequencies around the O2 A band of oxygen, is fairly constant 
around 10/m, i.e. of the order of 2.5 x 10^ times stronger than the absorption 
in O2 [24]. 

For the absorption coefficient in O2, we note that the 2.6% damping over 
one meter in the framework of Eq. (1) corresponds to ///q = 1 — 0.026 = 
exp(— 0.20942 x x 1 m) where we have explicitly accounted for the volume 
fraction of oxygen being 20.942% in ambient air. This translates into = 
0.126 m"^ and for PMMA we have ctPMMA = 10 m ^ while the corresponding 
dielectric functions are e'o^ = 1^ and Cpj^jviA = 1-5^. Choosing A = 1 /xm we 
find that the oxygen A band will be located right next to a band gap edge if 
we choose hg = 0.6262A, resembling the situation shown in Fig. 4. Assuming a 
Lorentzian absorption profile and using the method discussed above, we find 
that 702 — 235. We emphasize that the Lorentzian linewidth correspond to 
the complete A band and thus resembles a group of different and even sharper 
absorption lines, see inset of Fig. 5, for which the enhance factor varies. 

In a typical Beer-Lambert experiments measuring the O2 concentration in 
the atmosphere, a path length of around 100 cm gives a 2.6% attenuation 
of the transmitted signal. The enhanced absorption in the Bragg stack dis- 



9 



cussed above implies that one could scale the path length down by a factor 
of 235^ ~ 145 to about 7 mm while retaining the same absorption signal. 
We emphasize that this sensitivity enhancement factor (effective path length/ 
real length of the absorber) is even higher than that of a multi reflection cell, 
which is limited by reflection losses to about 100! Very high enhancement fac- 
tors of 5 orders-of-magnitude or more compared to a single pass setup may be 
achieved using cavity-enhanced or cavity ringdown laser spectroscopy. How- 
ever, such macrooptical setups can hardly be condensed into a microoptical 
device, and are generally delicate instruments mainly operated in laboratory 
environments. Using pairs of highly reflecting mirrors with reflectivities ex- 
ceeding 99.999%, these setups are very sensitive to contamination e.g. by 
dust. In contrast the slow-light concept can be considered as a manifold of 
distributed mirrors. Hopefully, contamination shall not be as detrimental as 
with the macroscopic cavity approaches. In terms of bandwidth, the slow- 
light phenomenon shares physics with the various cavity-based approaches: 
strong enhancement of the absorption is obtained at the price of a narrow 
band-width! 



6 Conclusion 

Slow-light phenomena in photonic crystals has recently been proposed as a 
mechanism to allow for reduced path length in absorption measurements in 
microsystems. In this paper we have considered the effect an absorption line 
in a gas infiltrating a uniformly absorbing Bragg stack. We have shown that 
enhanced absorption is possible and as an example we have considered a Bragg 
stack consisting of PMMA infiltrated by O2. In that case we have found the 
absorption enhancement to exceed a factor of 10^ potentially allowing for O2 
detection with sub millimeter optical path lengths. For future directions we 
emphasize studies of the influence of the interfaces (scattering loss), struc- 
tural imperfections, and beam divergence in more detail and another natural 
extension would be planar two-dimensional photonic crystals membranes. 
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